Theory of Distinct Crystal Structures of Polymerized Fullerides ACeo? A=K, Rb, Cs: 

the Specific Role of Alkalis 
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The polymer phases of AC60 form distinct crystal structures characterized by the mutual orienta- 
tions of the {Cq )„ chains. We show that the direct electric quadrupole interaction between chains 
always favors the orthorhombic structure Pmnn with alternating chain orientations. However the 
specific quadrupolar polarizability of the alkali metal ions leads to an indirect interchain coupling 
which favors the monoclinic structure 12 /m with equal chain orientations. The competition between 
direct and indirect interactions explains the structural difference between KCeo and RbCgo, CsCeo- 
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Alkali metal doped Ceo (A^Cgo), A=K, Rb, Cs, form 
stable crystalline phases (fullerides) with a broad range 
of physical and chemical properties comprising supercon- 
ductors and polymer phases. For a review, see jjJU. 
In particular the x = 1 compounds Q exhibit plas- 
tic crystalline phases with cubic rocksalt structure at 
high temperature (T > 350 K) and polymeric phases 
of reduced symmetry at lower T. In the latter the 
Ceo molecules are linked through a [2+2] cycloaddition 
||, a mechanism originally proposed for photo-induced 
polymerization in pristine Ceo- From X-ray powder 
diffraction it was concluded that the crystal structure 
of both KCeo and RbCeo was orthorhombic (space group 
Pmnn). Polymerization occurs along the orthorhombic 
a axis (the former cubic [110] direction), where the ori- 
entation of the polymer chain is characterized by the 
angle \i of the planes of cycloaddition with the c axis. 
In the Pmnn structure (Fig. la), these orientations are 
alternatively \x and —fx, \u\ w 45° ± 5°. Notwithstand- 
ing this apparent structural similarity, the electronic and 
magnetic properties of KCeo on one hand and RbCeo, 
CsCgo on the other hand, were found to be very differ- 
ent p. ESR and optical conductivity data ||||] show 
that RbCeo and CsCeo exhibit a transition from a quasi- 
one dimensional metal to an insulating magnetic state 
near 50 K, while KCeo stays metallic and nonmagnetic 
at low T. NMR spectra also showed marked differences 
between KCeo and Rb-, Cs- polymers f|. The contradic- 
tion between similar crystalline structures and different 
electromagnetic properties was resolved by single crystal 
X-ray diffraction and diffuse scattering studies JlCj] . In- 
deed the polymer phases of KCeo and RbCeo are differ- 
ent. While the space group Pmnn is confirmed for KCeo, 
it is found that RbCeo is body centered monoclinic, with 
space group I2/m. In the latter structure, the polymer 
chains have the same orientation u (Fig. lb). Electronic 
band structure calculations for the 12/ m structure have 
shown the importance of transverse interchain coupling 
for RbCeo fn] ■ Recently [O] high-resolution synchrotron 
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FIG. 1. Crystal structures projected onto the crystallo- 
graphic (b,c) plane: (a) Pmnn, (b) 72/m. The thick bars 
represent the projection of the cycloaddition planes. Poly- 
merization occurs along a. The alkalis located in (b, c) planes 
and at ±o/2 are denoted by full (+) and empty (— ) circles. 



powder diffraction results demonstrated that CsCeo has 
the same structure as RbCeo- 

In this paper we study the mechanism which leads to 
the distinct polymer phases and we demonstrate the ac- 
tive role of the distinctive quadrupolar electronic polar- 
izability of the alkali ions (cations). We start from the 
high temperature orientationally disordered cubic phase 
(space group Fm3m) and describe the formation of the 
polymer phases as a scenario with several steps: i) the 
charge transfer of one electron from the alkali atom to 
the Ceo molecule leads to an occupation of the lowest 
unoccupied molecular orbital levels which are of t\ u sym- 
metry. Thereby the crystal field of the Ce~ ion acquires 
an electronic component [fl3| which favors a same ori- 
entation of neighboring molecules along [110] such that 
the stereospecific cycloaddition occurs, ii) The cycload- 
dition between neighboring molecules then acts as a neg- 
ative internal stress (chemical pressure) along [110]. iii) 
The mutual orientation of neighboring C^~ chains, which 
distinguishes between the orthorhombic and monoclinic 
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structures, should depend on the intercalated alkali ions. 

We first show that the orthorhombic lattice is the re- 
sult of the cycloaddition. Using concepts of the theory of 
elasticity |Oj , we find that the cubic crystal is deformed 
into an orthorhombic one (point group D 2 h)- Taking 
the cubic [110], [110] and [001] as new x, y and z axes 
(orthorhombic a, b and c) respectively, we find the defor- 
mations 
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dc 44 
fzz = -Kcvz/d 



[cii(cn + C12 + 2c 44 ) - 2c\ 2 } 7 
[cn(cu + C12 - 2c 44 ) - 2cf 2 ], 



(la) 

(lb) 
(1c) 



where are the cubic elastic constants, d — Cn(cn + 
c 12 ) — 2cf 2 and where K < is the uniaxial stress. Obvi- 
ously e xx < 0, c yy > and e zz > which corresponds to 
contraction along a and elongations along b and c. We 
next investigate the orientation of the pol ym er chains. 
Following the experimental breakthrough |10|, prelimi- 
nary calculations of interchain energies have been per- 
formed |]l5| . The interaction energy is highly sensitive to 
the lattice constants and a plausible scenario is that Ceo- 
C6o interchain distances impose different relative chain 
orientations. However the last argument is at variance 
with single crystal X-ray diffraction results of pressure 
polymerized Cgo p6| . There the space group is Pmnn, 
i.e. isostructural with polymerized KCgo- However the 
orthorhombic cell volume of polymerized Cgo and the dis- 
tance between a corner and the center of the cell, 1326 A 3 
and 9.956 A, respectively, are closer to the corresponding 
values in RbC 60 , 1314.9 A 3 and 9.852 A, than to those of 
KC 60 , 1298 A 3 and 9.836 A. We conclude that the alkalis 
must play a more specific role in triggering the struc- 
tural difference between KCeo and RbCgo, CsCgo- From 
the theory of bilinear rotation-translation (RT) coupling 
between molecular rotations and lattice displacements of 
the counterions which plays an essential role in determin- 
ing the elastic properties of ionic molecular crystals |L7j , 
one finds an effective rotation-rotation (RR) interaction 
that competes with the direct intermolecular RR inter- 
action. However the lattice mediated RR interaction is 
found to be independent of the alkali mass and hence no 
distinction on basis of the different masses is possible. A 
further distinctive property of the alkali metal ions is the 
dipolar electronic polarizability, with values 0.9, 1.7 and 
2.5 A 3 for K+, Rb+ and Cs+, respectively We ex- 
pect that also the quadrupolar polarizability is larger for 
the heavy alkali metal ions than for K + . In the following 
we will show that the electric quadrupole interaction be- 
tween C^g chains and the alkali ions leads to an effective 
orientational interaction between the C^~ chains. 

We will study first the electric multipole interaction be- 
tween polymer chains. Polymerization reduces the sym- 
metry, the C^q chain is composed of units with D 2 h sym- 



metry. We have used a tight-binding model to study 
the electronic charge distribution on the Cg units in 
the chain [ jl9| . The charge is mainly concentrated in 
the equatorial region of Cqq, in agreement with recent 
NMR results p0[ ]. We find that only even I multipoles 
are allowed, in particular each C^~ unit has an electric 
quadrupole. In the following we adopt a simple model of 
charge distribution. Using the labeling of C atoms of [^0| , 
we locate a charge of —0.15 (units |e| = 1) on each bond 
C15-C16. These charges are fixed at a distance d = 3.52 
A from the center of the Cg ball. Such a charge dis- 
tribution is sufficient to obtain a quadrupole (see also 
Eq. (|^)), the accompanying monopole is irrelevant. The 
chains are taken as rigid units with sole degree of free- 
dom the rotation angle fi about the axis a. Hence we 
can treat the three dimensional crystal in the polymer 
phase as a two dimensional problem in the (be) plane. 
Here we consider one chain per unit cell with basis vec- 
tors f x = (c/2) + (6/2) andr 2 = (c/2)-(6/2). The chains 
are labeled by a two dimensional array ft — (nx,ri2) 
and lattice vectors are given by X(n) = nxfx + nif-x-, 
where n\, n% are integers. The Coulomb interaction be- 
tween chains depends on their mutual orientation. We 
introduce symmetry adapted rotator functions (SAF's) 
Si(n) = sin(//i(r?)). Symmetry of the chain for /i — > /i + 7r 
implies I — 2, 4, ... . These functions are ungerade in /j, 
and S*2 is maximum for /i = 45°. Expanding the Coulomb 
interaction in terms of SAF's we consider in lowest order 
of I the quadrupole-quadrupole interaction 



H S s = \ y^(J a (n,n') + J b (H,H')) 5 2 (n) S 2 (n). 



(2) 



More involved molecular charge distributions would lead 
to higher (I > 4) multipoles |I3] ], but their interactions 
are negligible. We observe that nearest neighbor chains 
are shifted by a translation a/2 along the a axis with re- 
spect to each other. For a C^ molecule at n as origin, J a 
takes into account the interaction with the two molecules 
at ±a/2 on the neighboring four chains X(n') — ±rj 
and X(n') = ±r 2 , <h describes the interaction with one 
molecule on the chains X(ft') = ±(fi — f%). In two di- 
mensional Fourier space we obtain 

J(q) = 8J a cos (m cos (^pj + 2J b cos( % 6), (3b) 



where q y and q z are components along the original or- 
thorhombic axes b and c. With the simple model of 
charge distribution for the C^~ units and the orthorhom- 
bic lattice constants of Table I, we have calculated the 
quadrupolar interaction energies J a and ,J\, quoted in 
Table I. From Eq. (3b) we find that J(q) is negative 
and maximum in absolute value at the boundary of the 
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TABLE I. Lattice constants a c refer to cubic, a , b and 
c to simple orthorhombic lattice cells, length in units A 
(Angstrom). Interaction energies J a , Jb, A5, A c , qa measured 
in units K (Kelvin). Quadrupolar radii (La in units A. Ratios 
of qa 's are equal to the corresponding ratios of I/^a's. 



H sS = }X\b(n,n') + X c (ft,n')) S 2 (n) s 2 (n'), 



(4) 
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Brillouin zone (BZ). For all three compounds we have 
|J(<fz)| > l^(Sr)li w here qz — (0, 2-k/c) is a BZ boundary 
vector and qr — (0, 0) is the BZ center. This result is in- 
dependent of the strength of the molecular quadrupoles 
but is a consequence of the orthorhombic lattice. The 
dominance of J(qz) leads to a condensation of S 2 (q) for 
<f = Qz- (5*2 (<lz)) = "q/y/W. Here 77 is the order parame- 
ter amplitude and N is the number of chains (i.e. lattice 
points in the (be) plane). Condensation at qz implies that 
the chains in a same basal plane (ab) of the orthorhombic 
lattice all have the same orientation, but the orientation 
alternates in neighboring planes at distance c/2. This is 
the "antiferrorotational" structure Pmnn, Fig. la. We 
find that the quadrupolar interaction between Cgo chains 
leads to Pmnn for KCgo, RbCeo, CsCeo, irrespective of 
the different orthorhombic lattice constants. For RbCgo 
and CsC 60 the experimental result is 12 /m pX|Jl^ |. We 
now include the role of the alkali ions with their distinc- 
tive polarizability. It is known from work on the ammo- 
nium halides H] that the indirect interaction of two NH4" 
tetrahedra via the polarizable halide ions plays an essen- 
tial role in determining the various crystalline phases of 
the ammonium halides NH4X, X=C1, Br, I p2[ . However 
in the present problem, by symmetry the dipolar polar- 
izability of the alkali metal ions is irrelevant and we have 
to resort to the quadrupolar polarizability. Since the C^~ 
units in a polymer chain are rigidly linked in a same ori- 
entation, the C^q chains produce coherent electric field 
gradients which induce an anisotropic (quadrupolar) de- 
formation of the electron shell of the alkalis. We model 
the corresponding charge distribution of each alkali by 
a symmetric linear dumbbell centered on lines along a 
directions. On a same line, these induced dumbbells are 
parallel with their axis perpendicular to a and a same ori- 
entation angle v with the c axis. We consider chains of 
alkali dumbbells, where the rigid chain aspect is not im- 
posed by intrachain interactions but by the surrounding 
C^g chains. The orientational motion of an alkali dumb- 
bell is characterized by the SAF's s 2 (n) = sin(2^(n)), 
where ft is again a two dimensional array labeling the 
chains. The quadrupole-quadrupole interaction between 
the C^g chains and the surrounding alkalis is then given 



by 



where Xb(n,n') accounts for the two alkalis at ±a/2 on 
the chains X(n') = ±(r*i — ^)/2 and where X c (n,n') 
describe the interaction with one alkali on the chains 
X(n') = ±(r*i + r*2)/2. In Fourier space we have 

H sS =Y / X(q)S 2 (q)s 2 (-q) 7 (5a) 

X(q) = 4A b cos + 2X C cos (^) . (5b) 

We take dumbbells with charges qA at distances ±g?a 
from the center. Here A refers to K, Rb, or Cs. The 
numerical values of d,A (Table I) are the average radii 
of valence electron d shells calculated with atomic wave 



functions 3c? 



3/2, 4o! 3 / 2 , 



5o! 3/2 for K+, Rb+ and Cs+ re- 



spectively. We consider d-shells because they can support 
an electric quadrupole moment. We observe that the val- 
ues for Cs + and Rb + are close to each other but differ 
from K + . The alkali ions are isoelectronic with the rare 
gas atoms Ar, Kr, Xe. There the role of excited d states 
has been found to be important in the explanation of the 
face centered cubic structure |E3[. With the same value 
qA = 0. 15 |e I for the three cases, we have calculated the 
interaction energies Ab and A c quoted in Table I. The 
quantity |A(<t*)| is maximum for q = qr in contradistinc- 
tion with |J(g)|. The intra ionic restoring forces which 
oppose the deformation of the electron shells of the alkalis 
are described by a sum of single particle energy terms 



Hss = 9A^2,s%(n) = g A ^s 2 (q) s 2 (-q), 



(6) 



with gA > 0. The self energy gA is inversely propor- 
tional to the quadrupolar electronic polarizability and 
hence gc s < gsb < gn (sec Table I). These concepts are 
inspired from the shell model of lattice dynamics 
The direct interchain coupling of alkali quadrupoles is 
numerically found to be small and will be neglected. We 
consider the total interaction H = Hss + H s s + H ss . The 
induced alkali quadrupoles follow adiabatically the mo- 
tion of the C^q chains. For a given configuration {S 2 (q}} 
of the latter, we minimize H with respect to s 2 (q) and 
find 



s 2(q) 



2 g A 



S 2 (q), 



Substituting this result for s 2 in H we get 
H = \Y,J(q)S 2 (q)S 2 (-q) 1 



where 



J(q) = J(q) + C(q), 

l A 2 (g) 
2 g A 



C(q) 



(7) 

(8a) 

(8b) 
(8c) 
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FIG. 2. Interchain energy J(q), units K. 



The coupling to the alkalis leads to an effective orienta- 
tional interaction C(q) between C^" chains. This inter- 
action is attractive and since |C(q*)| with \{q) is max- 
imum at q = qr, it favors a "ferrorotational" struc- 
ture (space group J2/m). On the other hand the di- 
rect quadrupolar interaction between C^~ chains with 
|J(<fz)| > |J((fr)| favors the antiferrorotational struc- 
ture. The system of interacting polymer chains described 
by Eq. (^) chooses the lowest energy structure which is 
Pmnn or I2/m depending on whether | J(gz)\ > \ J(Qv)\, 
i.e. J a > Xb^c/dA, or the opposite holds respectively. Us- 
ing the numerical values from Table I, we have plotted 
J(q) for q=(q y = 0, q z ) in Fig. 2. For KC 60 § with 
the small polarizability of the K + ion, the direct inter- 
chain interaction Jiqz) dominates and leads to Pmnn 
while for RbC 60 Q and even more for CsC 60 JlJ] the 
alkali mediated interaction C(qp) dominates and leads to 
12 /m. A condensation of S-2,(q) at the BZ center leads, 
via coupling to the center of mass displacements of the 
alkali ions (bilinear RT-coupling) to e yz shear modes and 
hence to a deviation of the (b, c) angle a from 90° in 
the 12 /m structure. However this RT-coupling is not the 
driving process of the monoclinic structure. Indeed, the 
deviations of a are very small jnj[^ . 

In conclusion we have shown that the cycloaddition be- 
tween C^q units leads to an orthorhombic lattice struc- 
ture. The concomitant symmetry reduction produces a 
quadrupolar electric charge distribution on the C^q units. 
The interaction between C^" chains has two competing 
components: a direct quadrupole-quadrupolc interaction 
J(q) and an indirect one C{q). The latter is mediated 
by the induced quadrupoles on the polarized alkali metal 
ions. The direct quadrupolar interaction drives the an- 
tiferrorotational structure Pmnn while the indirect one 
yields the ferrorotational structure 12 /m. In RbCgo and 
CsCgo with larger electronic polarizability and quadrupo- 
lar radius oIa of the alkalis, XbX c /gA > Ja, the ferrorota- 
tional structure is realized while in KCgo, J a > Xb\ c /gA, 



the antiferrorotational structure is realized. Within the 
present theory, the alkalis play a specific role beyond the 
function of lattice spacers. The study of alkali specific 
effects is a problem of broad interest and likely to be rel- 
evant for the understanding of superconducting fullerides 
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